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Introduction 

Un treillis est un ensemble de poutres droites (éléments) reliées entre elles par des rotules (nœuds). 
Les liaisons extérieures sont des rotules et des appuis simples. Les charges sont des forces portées par 
les rotules, des gradients thermiques et des déplacements d’appui. La force intérieure dans une 
section droite se réduit à l’effort normal. 

Le treillis est plan si : 

Le plan {0;x,y} est un plan de symétrie pour toutes les sections droites. 

Les forces appliquées sont situées dans le plan {0]x, y}. 

On suppose que les déplacements sont petits. 


1 Matrices élémentaires 

Soit (i —> j) un élément de treillis plan de section droite constante (figure H]). 




y 


► X 



Figure 1 - Elément i —>• j 

L est la longueur de l’élément et A l’aire de sa section droite. 

(xi , iji) et ( Xj , 'ijj) sont les coordonnées des nœuds de l’élément. 
Le vecteur unitaire n porté par l’axe de la poutre est défini par : 


1 txj - ( X , - Xl f + ( Vj - 


n x 

n v) L l V 3 ~ Vi 


( 1 . 1 ) 


où 6 est l’angle que fait n avec l’axe x. 

E et a sont respectivement le module de Young et le coefficient de dilatation du matériau. 


L’élément est soumis à un effort normal N (positif : traction, négatif : compression) et à une variation 
de température AT constante. 
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Méthode des éléments finis 


( Ui,vi ) et ( Uj,Vj) sont les déplacements nodaux (figure EJ). 





► x 


Figure 2 - Déplacements élémentaires 

Les efforts aux extrémités de l’élément sont (figure [3]) : 

— N n en i , N h en j 


( 1 . 2 ) 



N n 


Différentions la relation 

Il vient : 


Figure 3 - Efforts élémentaires 
L 2 = ( Xj - Xi) 2 + (yj - yf) 2 


2 L dL = 2 (xj - Xi) ( dxj - dxi) + 2 {yj - yf) (dyj - dyf) 

d’où l’expression de l’allongement unitaire suivant n : 

dL 1 f (xj -xf) {yj - yf) , 

En. = — = — I — — - {dxj - dxi) H- J — - {dyj - dy i/ 


soit : 


L L V L 


1 


E n = ~{ n x {uj - m) + n y (vj - Vj) ) 


(1.3) 

(1.4) 

(1.5) 

( 1 . 6 ) 


Cet allongement unitaire est dû à l’effort normal (loi de Hooke) et à la variation de température : 

N 


e n - — +aAT 

L’effort normal s’écrit en fonction des déplacements nodaux : 


N = EA (s n — a AT) 
_ EA 


( n x {uj — uf) + n y (vj — Vi) ) — EA a AT 


(1.7) 


( 1 . 8 ) 


soit : 


EA 




n x . riy n x n y J < 


Ui 

Vi 

Ui 


- EAaAT 


(1.9) 









Treillis plans à nœuds articulés 
On en déduit : 


{/nod} = [fc] M - {Ah} 


avec : 


' -n x 


' u t \ 


II 

's' 

Vi \ 

n x 

u 0 ! 

. n y > 


. v j ) 


EA 

~L~ 


’ ~nx 



‘ n l 

flx fly 

-ni 

Tlx Tly 

n x 

. n y > 

> [~n x 

-, EA 

Tiy Tlx Tlyj — 

sym. 

n 2 

Uy 

Tlx Tl y 

nî 

— T? 

n y 

Tlx Tly 

n? 

n y J 


~n x 

{fth} = EAaAT {~™y 

,I 'X 

n,, 


1 y 

{/nod} est le vecteur force nodal (N). 

{u\ est le vecteur déplacement élémentaire (m). 

[ k ] est la matrice de rigidité (N/m). 

{Ah} est le vecteur force équivalent au gradient thermique (N). 


Remarque 1 : la matrice de rigidité peut se mettre sous la forme : 
[k] = 


M -[ k} 

-m [ à ] 


m EA 

avec [k\ = —— 

1 j 


'°x ' ü x '°y 

n x n y n y j 


Remarque 2 : la matrice de rigidité est égale à : 

[*] = 


n x 

0' 





n y 

0 

EA 

' 1 

-1' 


n x n y 0 0 

0 

n x 

~L 

-1 

1 


0 0 n x n y 

0 

n y _ 






ou : 


[k) = 


n x 

-% 

0 

0 


' 1 

0 

-1 

0' 


n x 

n v 

0 

o' 

n y 

n x 

0 

0 

EA 

0 

0 

0 

0 


~n v 

flX 

0 

0 

0 

0 

n x 

—n y 

~L~ 

-1 

0 

1 

0 


0 

0 

Tlx 

n y 

0 

0 

ny 

n x _ 


0 

0 

0 

0 


0 

0 

-Tly 

n x _ 


(1.10a) 


(1.10b) 


(1.10c) 


(l.lOd) 


( 1 . 11 ) 


( 1 . 12 ) 


(1.13) 


Remarque 3 : l’énergie de déformation est égale à (à un coefficient près indépendant des déplacements 
et de leurs dérivées) : 


Edeî — 2 EA s^ L — EA s n a AT L 

= \ M T [k] M - M T {/th} 


(1.14) 


Le travail des forces extérieures se réduit au travail des forces nodales : 


Wext = M T {/nod} 


( 1 . 15 ) 
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Méthode des éléments finis 


L’énergie potentielle est égale à : 


■Epot — Edef ILext 


(1.16) 


La matrice de rigidité est la matrice hessienne (ou matrice de Hess) de l’énergie de déformation par 
rapport aux déplacements nodaux (programme ltre_matl) : 


_ <9 2 Edef 

y dui duj 


(= kji ) 


(1.17) 


Le vecteur des efforts aux nœuds est le gradient de l’énergie de déformation par rapport aux dépla¬ 
cements nodaux : 


fnod,i 


dE dei 

ÔUi 


(1.18) 


2 Exemple 1 : treillis soumis à une force nodale 

2.1 Énoncé 


Le treillis plan à noeuds articulés représenté sur la figure [4] est com¬ 
posé de trois poutres de même nature et de même section droite. 

Soient E le module de Young du matériau et A l’aire des sections 
droites. 

Le noeud 1 est articulé et le nœud 3 repose sur un appui simple dont 
la normale est horizontale. 

Le noeud 2 porte une charge de composantes (0,P). 

Application numérique : on donne : 



A = 100 mm 2 ,L = 0.2m,E = 200000 MPa , P = -10000 N Figure 4 Exemple 1 


2.2 Partition des degrés de liberté 

Effectuons une partition des degrés de liberté en déplacements connus et inconnus (m m) ■■ 


{u L } = |v2 J 



d’où {U} 



' V,2 ' 
V2 
V3 
Ul 
Vi 
.«3, 


On en déduit la localisation des degrés de liberté 


dans les matrices globales : 


'm -> 0' 
ui —y 0 
U 2 ~> 1 
V'2 —> 2 

«3 —> 0 

. ^3 t 3 j 


{DDL} 
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2.3 Étude élémentaire 

- coordonnées nodales : 


élément 1 —> 2 : 


nœud 

X 

y 

1 

0 

L 

2 

L 

0 

3 

0 

- L 


caractéristiques : L \[2 , A , E , n x = — y= , n y = — —^= 


fui-xn 


' 1 -1-1 1 ' 

v\ —> 0 

r, 1 EA 

-11 1-1 

U2 1 

’ [h - 2] -2V2 L 

-11 1-1 

U 2 ^2) 


. 1 -1-1 1 . 


élément 3 —> 1 : 


caractéristiques 

: 2 L 

A, E 

; Tlx — 

0 , 

n y = 

1 




-)• 0' 




0 

0 

0 

0 ' 

i}=< 

V 3 

3 


[^ 3 - 1 ] 

EA 

0 

1 

0 

-1 

U\ 

-> 0 

> , 

~ 2T 

0 

0 

0 

0 


Ai 

-)• 0 

/ 




0 

-1 

0 

1 _ 


élément 3 —» 2 : 


caractéristiques : L \]~2 , A , E , n x = —= , = —= 

v 2 y 2 


{<M3_2} — < 



) [^ 3 - 2 ] = 


EA 
2y/2 L 


1 1-1 - 1 ' 
1 1-1-1 
-1-11 1 

-1-11 1 


2.4 Assemblage et calcul des déplacements inconnus 

Les déplacements inconnus sont les solutions de l’équation [Kll]{Ul] = {Éiod,i} : 


EA 


2y/2 L 


'2 0 -1 
0 2-1 
-1 -1 l + y/2 


U2 

^2 ) = { P 

V 3 


d’où (programme |exemple_l] ) : 

PL 


PL 


u 2 = = -0.050 mm , v 2 = (1 + 2\/2) = -0.191mm 

Z JbA Z h/A. 

PL 

v 3 = -^-r = —0.100 mm 


EA 

2.5 Efforts normaux dans les éléments 

Ils sont calculés à l’aide de la formule 


N - = l|(^-^ = -^ =7071N 

EA P 

N 3-1 = jï {~v 3 ) =-j = 5000N 


A 3—2 = 


EA ( 1 


1 


—7= U2 + —7= ( V 2 - V 3 ) = 


Ly/2 W2 y/2 


Psj2 


= -7071N 
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2.6 Actions de liaison 

Les actions de liaisons sont calculées à partir des efforts normaux (équation 11 .21) : 
nœud 1 : 

Fi = — Ai _2 ni -2 + A 3 _i n 3 _i d’où 


Fix 


—\= Ai_ 2 = — = -5000 N 
y/2 2 


Fl y ~ 


1 

71 


Ai_2 + A 3 _i = —P = 10000 N 


nœud 3 : 


F 3 = -JV 3 _i n 3 _i - A 3 _ 2 n 3 - 2 d’où 
F 3x = —— A 3 _2 = — — = 5000 N 


Remarque : l’équilibre de la structure est vérifié : 


Fix + F2x + F 3æ — 0 , F\y + F2y + F 3y — 0 , —2 LF\ X — LF2 X + LF2 y — 0 


3 Exemple 2 : treillis soumis à une force nodale 


3.1 Énoncé 


Le treillis plan représenté sur la figure [5] est composé de trois 
poutres de même section. 

Soient E le module de Young du matériau et A l’aire des sec¬ 
tions droites. 

Le nœud 1 est articulé et le nœud 2 repose sur un appui simple 
dont la normale est horizontale. 

Le nœud 3 porte une charge d’intensité (F, 3 P, 0). 

Application numérique : on donne : 



P 

3P 


Figure 5 - Exemple 2 


A = 100 cm 2 , E = 200000 MPa , L = 0.7 m , P = -120 kN 


3.2 Partition des degrés de liberté 

Effectuons une partition des degrés de liberté en déplacements connus et inconnus (ffUEl): 
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On en déduit la localisation des degrés de liberté dans les matrices globales : 


{DDL} 


' u\ — y 0 ' 
v\ — > 0 


u 2 —> o 

V2 —> 1 




U 3 —> 2 

A3 3 _ 


3.3 Étude élémentaire 

- coordonnées nodales : 


élément 1 —>• 2 : 


nœud 

X 

y 

1 

0 

0 

2 

0 

L 

3 

L 

0 


caractéristiques : L , A , E , n x = 0 , n v = 1 


' U\ 

0' 



0 

0 

0 

0 ' 

V\ 

0 


ri 1 EA 

0 

1 

0 

-1 

U2 

-A 0 

> 

[h-2] - -j- 

0 

0 

0 

0 

K V 2 

-)• 1 

/ 



0 

-1 

0 

1 


élément 1 —> 3 : 


caractéristiques : L , A , E , n x = 1 , n y = 0 

u\ —> Ol 1 

ii._i.nl K A n 

{ddli-3} = < 


'ui -f 0Ï 

' 1 

0 

-1 

0' 

V\ —> 0 1 r. , EA 

0 

0 

0 

0 

«3^2 • |fcl - 31 = ~L 

-1 

0 

1 

0 

, v 3 -> 3 J 

_ 0 

0 

0 

0 


élément 2 —y 3 : 


caractéristiques : L y/2 , A , E , n x = —= , n y = —= 

V 2 y 2 


{(M 2 - 3 } — < 


fn 2 ^o] 


' 1 -1-1 1 ' 

V2 — > 1 

n ! EA 

-11 1-1 

u-s -> 2 

’ [ 2 ~ 3] ~2y/2 L 

-11 1-1 

^ 17 3 — ► 3 J 


1-1-11 


3.4 Assemblage et calcul des déplacements inconnus 

Les déplacements inconnus sont les solutions de l’équation [Kll\{Ul} = {Éiod,L} : 


EA 
2 y/2 L 


l + 2y/2 1 -1' 

1 1 + 2 y/2 -1 

-1 -1 1 


V2 


d’où (programme exemple_2) : 


3 PL 

V 2 = — . = —0.126 mm 


EA 


APL 

«3 = = —0.168 mm 

EA 


(7 + 6^2 )PL 

vs = -——-= —U.ooU mm 

EA 



























10 
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3.5 Efforts normaux dans les éléments 

Ils sont calculés à l’aide de la formule (|1.8D : 

EA 


N\—2 = —— V 2 = 3 P = —360 kN 

1j 

EA 

1Vi_3 = ——— u 3 = 4 P = —480 kN 

1j 


N2-3 = 


EA f 1 


1 


Ly/2 W2 U3 s/2 

3.6 Actions de liaison 

Les actions de liaisons sont calculées à partir des efforts normaux (équation 11.21) : 

- nœud 1 : 

Fi = -lVi_2 ni _ 2 - Ai _ 3 ni _ 3 d’où 
F\ x = —lVi_3 = —4-P = 480 kN , F\y = —fVi_2 = —3 P = 360 kN 

nœud 2 : 

F 2 = fVi -2 ni -2 ~ -ZV 2-3 n 2—3 d’où 
F 2x = — N 2-3 = 3 P = -360 kN 

s /2 

Remarque : l’équilibre de la structure est vérifié : 

Fix + Fïx + F 3 X = 0 , F\ y + F2 y + F^y = 0 , LF 3y — LF 2 . r = 0 


( v 3 - v 2 ) ^ = —3\/2 P = 509 kN 


4 Exemple 3 : treillis soumis à une variation de température 

4.1 Énoncé 

Le treillis plan à noeuds articulés représenté sur la figure [6] est composé de trois poutres de même 
matériau et de même section droite. 



Figure 6 - Exemple 3 

Soient E et a respectivement le module de Young et le coefficient de dilatation du matériau. 


Soit A l’aire des sections droites. 
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Les nœuds 1, 2 et 4 sont liés à l’extérieur par une rotule. 

Premier cas de charge : la structure est soumise à une variation de température AT. 
Deuxième cas de charge : la poutre (2 — 3) est soumise à une variation de température AT. 
Application numérique : on donne : 

A = 100 mm 2 , L = 0.1 m , E = 200000 MPa , a = 10~ 5 KT 1 , AT = 100 K 


4.2 Partition des degrés de liberté 

Effectuons une partition des degrés de liberté en déplacements connus et inconnus (cuni): 



{Us} 


'ui' 

Vi 

U 2 

V2 

U4 

.VA, 


d’où {U} 



'm' 

V3 

Ul 

VI 

U2 

V2 

U 4 

.VA, 


On en déduit la localisation des degrés de liberté dans les matrices globales : 


{DDL} 


' u\ -> 0 
va —y 0 
U 2 ~> 0 
V 2 ~> 0 
U3 t 1 




V 3 2 
ua —> 0 
. va t 0 _ 


4.3 Étude élémentaire 

- coordonnées nodales : 


élément 1—^3: 


nœud 

X 

y 

1 

0 

0 

2 

L 

0 

3 

L 

L 

4 

2 L 

L 


caractéristiques : L \[2 , A , E , a , n x = —^= , n y = — y= 


{ddh- 3 } — < 



[h-s] = 


EA 
2 y/2 L 


1 1 - 1-1 

1 1 - 1-1 

-1-11 1 

-1-11 1 


{/th,i— 3 } — EAa AT < 



caractéristiques : L , A , E , a , n x = 0 , n y = 1 
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{(M 2 - 3 } = < 


112 -> 0 
v 2 — > 0 

u;i -ï 1 

1^3 2 


r; 1 EA 

[fc2— 3 ] = 


0 0 0 0 ' 

0 10-1 
0 0 0 0 

0-101 


{/th,2— 3 } = EAaAT < 


0 

-1 

0 

1 


élément 3 —)■ 4 : 


{cM 3 _ 4 } = < 


caractéristiques : L , A , E , a , n x = 1 , n y = 0 

U3 —> 1 r 1 0—10' 

w 3 2 I . EA 0 0 0 0 

u A 0 ^ 3_4 ^ — ^zr —1 0 1 0 

v A —> 0 0 0 0 0 


{/th, 3 — 4 } = EA a AT < 


f-lï 

0 
1 
0 




4.4 Premier cas de charge 

4.4.1 Assemblage et calcul des déplacements inconnus 

Les déplacements inconnus sont les solutions de l’équation [. Kll\{XJl } = {Pth,i} : 


EA 


1 + 2V2 


1 


U3[ =' EA ^T^l 2 2 


mm 


2^2 L [ 1 1 + 2 V 2 J \v 3 

d’où (programme |exemple_3 ) : 

u 3 = (y/2 — 2) L a AT = — 0.0586mm , v 3 = V2 La AT = 0.1414 

4.4.2 Efforts normaux dans les éléments 

Ils sont calculés à l’aide de la formule (11.811 : 

fVi _ 3 =-— ( —= « 3 H —— v 3 ] — EA a AT = (\/2 — 2) EAa AT = —11716N 

L y 2 \ \/2 v 2 J 

EA 

N 2-3 = -—v3- EAaAT = (y/2 - 1) EAaAT = 8284N 
JL 

FA 

1V 3 _ 4 = — (-It 3 ) — EAaAT = (1 - V~2) EAaAT = -8284 N 
JL 

4.4.3 Actions de liaison 

Les actions de liaisons sont calculées à partir des efforts normaux (équation 11,2|) : 
nœud 1 : 

F\ = — Ai _ 3 ni _ 3 d’où : 

F lx = —L Ai _ 3 = (y/2 - 1) EA a AT = 8284 N , F ly = F\ x = 8284 N 
v 2 

nœud 2 : 

F 2 = -N 2-3 T 2-3 d’où : 

F 2 X = 0 , F 21 J = —^ 2—3 = (1 — y/2) EAa AT = —8284N 

nœud 4 : 

F a = JV 3 _ 4 n 3 _4 d’où : 

E 4a; = A 3 _ 4 = (1 - V2) EA a AT = -8284 N , F Ay = 0 
Remarque : l’équilibre de la structure est vérifié : 

F\x + E 2 .X + Plx = 0 , F\y + + F A y = 0 , L F2y — L F Ax = 0 
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4.5 Deuxième cas de charge 

4.5.1 Assemblage et calcul des déplacements inconnus 

Les déplacements inconnus sont les solutions de l’équation [Kll\{Ul} = {Fl} '■ 


EA 


2 v/2 L 

d’où (programme exemple_3[ ) : 

1 - y /2 


1 + 2y/2 1 

1 1 + 2/2 


= EAot AT 


U 3 = 


La AT = —0.0207mm 


^ _ /?) 

V 3 = —-— L a AT = 0.0793 mm 


4.5.2 Efforts normaux dans les éléments 

Ils sont calculés à l’aide de la formule (11.81) : 

Nl - 3 = T7Ï (71 “ 3 + A ” 3 ) = ^ ^ “ Ar = 5858 N 

N 2-3 = —j— V 3 — EA a AT = - EA a AT = —4142 N 

L 2 

TT} A /Ô_ 1 

1V 3 _ 4 = — (-u 3 ) = -- EAaAT = 4142N 

L 2 


4.5.3 Actions de liaison 

Les actions de liaisons sont calculées à partir des efforts normaux (équation 11.21) : 


nœud 1 : 


1 


Fii ~~T 2 Ni ~ 3 ~ 2 


Fi = — fVi _3 ni _3 d’où : 
1 - y /2 


nœud 2 


F‘ 2 x = 0 , F 2 y = —N 2-3 = 


EAaAT = -4142 N , F ly = F lx = -4142 N 

F 2 = -N 2-3 n 2-3 d’où : 

y/ 2-1 


EAaAT = 4142 N 


nœud 4 : 


Fax — N3-4 — 


Fa = N 3—4 H 3-4 d’où : 

/Inl EAaAT = 4142 N , F iy = 0 


Remarque : l’équilibre de la structure est vérifié : 


Fia; + F 2 a; + Fax — 0 , Fly + F 2 y + Fly — 0 


L F 2y — L Fax — 0 
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5 Programmes Maple 

Les programmes suivant se trouvent dans le fichier treillis.txt. 

5.1 tre_mat 

# calculs élémentaires 
restart:with(linalg): 

# allongement unitaire 

eps:=(nx*(uj-ui)+ny*(vj-vi))/L; 


# énergie de déformation 

Edef:=EA*eps~2*L/2-eps*EA*alpha*DT*L; 

# matrice de rigidité 
k:=hessian(Edef,[ui,vi,uj,vj] ); 

# efforts nodaux 

fnod:=grad(Edef,[ui,vi,uj,vj] ); 

# remarque 

k:=j acobian(fnod,[ui,vi,uj,vj]); 

# vecteur dû au gradient thermique 
fth:=-jacobian(fnod,[DT]); 


5.2 exemple_l 

restart:with(linalg): 

# application numérique 

#L:=200 ;E:=200000 ; A :=Pi*30"2/4; P :=-10000 ; 

# matrice de rigidité 

KL :=matrix([[2,0,-1], [0.2.-1] ,[-1,-1,1+sqrt(2)]]): 
KL :=scalarmul(KL,E*A/2/sqrt(2)/L) ; 

# vecteur FL 

FL :=vector([0,P,0] ); 

# calcul des déplacements nodaux 
UL:=linsolve(KL,FL); 


#evalf (*/,) ; 


Treillis plans à nœuds articulés 
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5.3 exemple_2 

restart:with(linalg): 

# application numérique 

# L:=700 ;E:=200000 ; A :=10000 ; P :=-120e3; 

# matrice de rigidité 

KL :=matrix([[l+2*sqrt(2),1,-1],[1,l+2*sqrt(2),-l],[-l,-l,l]]): 
KL :=scalarmul(KL,E*A/2/sqrt(2)/L); 

# vecteur FL 

FL :=vector( [0,P,3*P]); 

# calcul des déplacements nodaux 
UL:=linsolve(KL,FL); 

#evalf(%); 


5.4 exemple_3 

restart:with(linalg): 

# application numérique 


#L:=100 ;E:=200000 ; A :=100 ; alpha :=le-5; DT :=100 ; 


# matrice de rigidité 

x:=l+2*sqrt(2): 

KL :=matrix([[x,1], [1,x]]): 

KL :=scalarmul(KL,E*A/2/sqrt(2)/L); 

# vecteurs F 

x:=E*A*alpha*DT/2: 

FLcasl:=vector([(sqrt(2)-2)*x,(sqrt(2)+2)*x]); 
FLcas2:=vector([0,E*A*alpha*DT]); 

# calcul des déplacements nodaux 

ULcasl:=linsolve(KL,FLcasl); 

ULcas2:=linsolve(KL,FLcas2); 

# map(evalf,ULcasl); 

# map(evalf,ULcas2); 
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